A two level finite difference scheme for one dimensional Pennes' bioheat equation

We develop a new two level finite difference scheme for the 1D Pennes bioheat equation. We further prove that the scheme is stable and convergent unconditionally. Numerical experiments for a skin-heating model are conducted. 2005 Elsevier Inc. All rights reserved. 0096-3003/$ see front matter 2005 Elsevier Inc. All rights reserved. doi:10.1016/j.amc.2005.01.052 * Corresponding author. E-mail addresses: [email protected] (J.J. Zhao), [email protected] (J. Zhang), nkang2@csr. uky.edu (N. Kang), [email protected] (F. Yang). URLs: http://www.umd.umich.edu/~xich (J.J. Zhao), http://www.cs.uky.edu/~jzhang (J. Zhang). 1 The research of this author was supported in part by the US National Science Foundation under grant CCR-0117602. 2 The research of this author was supported in part by the US National Science Foundation under grants CCR-9902022,CCR-9988165, CCR-0092532, and ACI-0202934, by the U.S. Department of Energy Office of Science under grant DE-FG02-02ER45961, by the Japanese Research Organization for Information Science and Technology, and by the University of Kentucky Research Committee. 3 The research of this author was supported in part by the US Department of Energy Office of Science under grant DE-FG02-0 ER45961. J.J. Zhao et al. / Appl. Math. Comput. 171 (2005) 320–331 321